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We present solutions to the equations of motion for bubble wall profiles in the 
minimal and a non minimal supersymmetric extension of the Standard model. 
We discuss the method of the numerical approach and present results for the two 
models (MSSM and NMSSM). 



1 Introduction 

For the emergence of a baryon asymmetry of the Universe the Sakharov condi- 
tions necessarily demand deviation from thermodynamical equilibrium. This 
condition is fulfilled in first order phase transitions. They take place via nuclc- 
ation of bubbles separating the symmetric from the broken phase. A first order 
phase transition might occur at temperatures around the electroweak scale. It 
turned out that in the Standard Model (SM) there is no phase transition at all 
for Higgs masses larger than 72 GeV. Baryon number generation at the elec- 
troweak scale therefore requires more complicated models with additional light 
scalar fields such as the MSSM or NMSSM. In the MSSM there is a window 
for electroweak baryogenesis and_an upper bound for the Higgs mass of about 
m H < 105 GeV with a light stojH of mass m t ~ H < m top . In the NMSSM the 
bound on the Higgs mass is even weaker □. 

Having established the existence of a first order phase transition one can 
start the actual calculation of the baryon asymmetry itself. There are several 
mechanisms described in the literature. All of them need the knowledge of the 
profile of the bubble wall during the phase transition. The kink ansatz in many 
situations is a good approximation but it might be interesting to have a more 
refined description and to determine which deviations occur in the presence 
of potentials depending on two or more Higgs fields and eventually on CP 
violating phases 0. Having the exact profile one can investigate the dynamics 
of expanding bubbles and calculate the baryon asymmetry. 

To determine the bubble wall profile beyond a simple ansatz we have to 
solve the equations of motion numerically. In the case of more than one scalar 



"Presented at "Strong and Electroweak Matter 98", Copenhagen, 2.-5.12.1998 
b Collaboration with S.J. Huber and M.G. Schmidt. 



1 



field this is a highly nontrivial task since simple methods like overshooting- 
undcrshooting fail. So one has to use methods which, beginning with an ansatz, 
converge to the actual solution. They are sometimes called "relaxation meth- 
ods". 

We first have to find the equations of motion. In field theory they can 
be derived via Euler Lagrange equations from the Lagrangian density which 
has the general form C = (D ll ^ i ) + (D^^ l ) + V($i,T) for several Higgs fields 
$i (plus, eventually, CP violating phases). Here Z? M is a covariant derivative 
and V denotes the effective potential. We show results of investigations of the 
MSSM and the NMSSM. 



1.1 Critical Bubble 

At the nucleation temperature the free energy becomes positive and the bub- 
bles start to grow until they fill up the entire space. The profile of the radial 
symmetric bubbles along the radius r is determined by the equations of motion 
for the critical bubble ("bounce"): 

Q= d 2 fa | idfa av , = 1 N 

dr 2 r dr dfa 

with boundary conditions = 0| r=oo and fa — 0| r=oo , where N — 2,3 for 
the MSSM and NMSSM respectively. 



1.2 Stationary Bubble or Domain Walls 

Constraining to a stationary wall with velocity v w at a late time t where the wall 
is already almost flat we are left with only one spatial dimension x = z — v w t 
perpendicular to the wall and have the equations 

with two boundary conditions, e. g. = 0|a; =oo and fa — 0| x= _oo. We 
will discuss now the method and show solutions calculated with it. 



2 Solutions 

2.1 The Method 

Since the overshooting-undershooting method fails we have to devise another 
method. We here minimize the functional of squared equations of motion. 
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Then, solving the eqs. of motion means finding field configurations for which 
the functional 

/+oo 
dx (E 2 (x) + E\{x) + . . . + E 2 N {x)) (3) 
-oo 

is zera, which is achieved by minimizing T . This method has already been 
used™ for the critical bubble of the MSSM. 

Applying the minimization method we have to solve a boundary value 
problem. Thus we have to make an ansatz for every function for which we 
want to find the time development which fulfills the boundary conditions. 

2.2 The Algorithm 

The numerical algorithm works in two steps: 

1. Find an ansatz as close as possible to the exact solution by low dimen- 
sional minimization of T with special ansatz configurations. 

2. High dimensional minimization of J- by discretizing every field function 
on a grid and minimization with respect to the function values. 

In the MSSM the kink ansatz is quite appropriate to get a good starting 
configuration for the minimization procedure of step 1. In the NMSSM the 
improved kink ansatz is not that good but nevertheless still appropriate. With 
only a few parameters Li and Xi this a low dimensional minimization procedure. 
In the MSSM we only have two (L\, L2), the offset x is negligible, in the 
NMSSM five parameters (£1/2/3, %i/2, x = 0). We use as ansatz configurations: 

</>f" fe = |(l + tanh(|- i = l...N (4) 

With N fields on a grid with M space points the second step is a N x M dimen- 
sional minimization, which must be performed with fast converging methods. 
The results here are obtained by N — 3 fields, M ~ 60 — 100 grid points and 
Powell's quadratically converging method. The derivatives of the differential 
equations are discretized with three and four point formulae, the integrals of 
T are performed with an. extended Simpson rule. The field configurations are 
interpolated by splines Em. For more details on the algorithm, see alsoO. The 
worst problem doing numerics is the existence of spurious minima. Besides 
real solutions to the equations there are fake minima due to the numerical 
representation and solutions due to the fact, that 8E 2 = is fulfilled not only 
for E = OJpjUt also for 8E — 0. One can perform checks to rate the minima 
found (seeH'Q). 
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Figure 1: Solutions for the MSSM. Left: Solution (solid) configuration and tan/i-ansatz 
(dashed). Only for unphysically small values of tua ~ O(10 GeV) the deviation from the 
straigth line (dashed) becomes appreciable. Wall widths arc around L f» 10/T C — 30/T c 
GeV" 1 . 



2.3 Applications 

We will show now applications of the described method. First, in figure 1 we 
present solutions of the bubble wall profile. The deviation from the straight 
line is small but nevertheless responsible for the actual amount of baryon asym- 
metry. Figure 2 shows the same for the three field case of the NMSSM. There 
we have a considerable deviation from the straight line and, additionally, a 
stronger deviation from the the extended kink ansatz. This demonstrates also 
the importance of a general solution method. Figure 2 shows also the path of 
the mechanical analogue of a rolling marble along the ridge of the potential. 
These results indicate the general behaviour of solutions in theories with more 
than one Higgs field. Now questions of metastability of wrong minima can be 
investigated with higher accuracy and better reliability. This is one step for a 
more precise calculation of the actual amount of the baryon asymmetry of the 
universe. 
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